Quasi-long-range order in trapped systems 
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We investigate the effects of a trapping space-dependent potential on the low-temperature quasi- 
long-range order phase of two-dimensional particle systems with a relevant U(l) symmetry, such 
as quantum atomic gases. We characterize the universal features of the trap-size dependence using 
scaling arguments. The resulting scenario is supported by numerical Monte Carlo simulations of a 
classical two-dimensional XY model with a space-dependent hopping parameter whose inhomogene- 
ity is analogous to that arising from the trapping potential in experiments of atomic gases. 
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Statistical systems are generally nonhomogeneous in 
nature, while homogeneous systems are often an ideal 
hmit of experimental conditions. Thus, in the study of 
critical phenomena, an important issue is how critical be- 
haviors develop in nonhoniogenous systems. Particularly 
interesting physical systems are interacting particles con- 
strained within a limited region of space by an external 
potential. This is a common feature of the experimental 
reahzations of the Bose-Einstein condensation (BEC) in 
diluted atomic vapors [l| and of optical lattices of cold 
atoms 01 , which have provided a great opportunity to in- 
vestigate the interplay between quantum and statistical 
behaviors in particle systems. 

In the BEC scenario, a macroscopic number of bosonic 
atoms accumulate in a single quantum state and are de- 
scribed by a condensate wave function, which naturally 
provides the complex order parameter ^{x) of the phase 
transition and its relevant U(l) symmetry. These global 
features characterize the XY universality class which de- 
scribes the universal critical behavior of a large class of 
systems, see, e.g., Ref. Q. The critical behavior aris- 
ing from the formation of the condensate in a trapped 
Bose gas has been investigated experimentally ob- 
serving an increasing correlation length compatible with 
the behavior expected at a continuous transition of ho- 
mogeneous systems belonging the three-dimensional (3D) 
XY universality class. Two-dimensional (2D) homoge- 
neous gases of bosonic particles do not show a real BEC 
with decreasing the temperature T . Neverthless, they 
are expected to experience a finite-T Kosterlitz-Thouless 
(KT) transition 5] separating the high-T phase with 
short-ranged correlations from a low-T phase character- 
ized by a quasi-long-range order (QLRO), where the one- 
body correlation function decays algebraically at large 
distance. Experimental evidences of such a transition 
in trapped Bose atomic gases have been provided in 
Refs. IsHSi. 

However, the inhomogeneity due to the trapping po- 



tential drastically changes the general features of the crit- 
ical behavior at the transition separating the high-T and 
low-T phases, and, in the case of 2D systems, of the 
QLRO phase. For example, correlation functions of the 
critical modes do not develop a diverging length scale in 
a trap. The critical behavior of the unconfined homo- 
geneous system could be observed around the middle of 
the trap only when the length scale f of the correlations 
is much smaller than the length scale induced by the 
trap size, and one looks at small-distance correlations rel- 
atively to . If ^ is large but not much smaller than the 
trap size, the critical behavior gets somehow distorted by 
the trap, although it gives rise to universal effects in the 
large trap-size limit, controlled by the universality class 
of the phase transition of the unconfined system [3, [13 ■ 
The understanding of the trap effects is necessary for an 
accurate determination of the critical parameters, see, 
e.g., Refs. dlHI. 

In this paper we consider 2D systems showing a low-T 
QLRO phase in their phase diagram, after a KT transi- 
tion. We investigate how the presence of the trap changes 
the main features of the QLRO of the homogeneous sys- 
tem and, therefore, how one may get evidence of the 
QLRO phase from the behavior of the system in the pres- 
ence of the trap. For this purpose, we resort to a scaling 
analysis which allows us to take into account the trap 
length scale when it becomes sufficiently large, exploit- 
ing the universality of the scaling behavior. 

The above considerations also apply to other physi- 
cally interesting models, such as the Bose-Hubbard (BH) 
model [l^ at its finite-T superfluid transition, whose 
Hamiltonian in the presence of confining potential reads 
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where the sum runs over the bonds (ij) of a d-dimensional 
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lattice, bi are bosonic operators, Ui = blbi is the parti- 
cle density, r the distance from the center of the trap. 
The trap size is defined as / = J^/^/w, see, e.g., [l4| . 
The BH model is of experimental relevance because it 
describes cold bosonic atoms trapped in a limited space 
region of optical lattices 13|. Other transitions in the 
XY universality class are the ''^He superfluid transition, 
insulator-superconductor transitions, like that of the at- 
tractive Hubbard model, etc.... 

In a standard scenario for a continuous transition, see, 
e.g., 0, the critical behavior of a d-dimensional system 
is characterized by two relevant parameters Ut and u^, 
which may be associated with T, i.e., ut ~ T/Tc — 1 
and the external field h coupled to the order parameter, 
with renormalization-group (RG) dimension yt — 1/v 
and Uh = {d + 2 — r])/2. The presence of a trap of size 
I generally induces a further length scale ^t, which must 
be taken into account to describe the critical correlations. 
Within the trap-size scaling (TSS) framework 0, the 
scaling law of the singular part of the free energy density 
around the center of the trap can be written as 



A standard representative model of the 2D XY univer- 
sality class is the classical square-lattice XY model. 



(2) 



where 9 is the trap exponent. At the critical point (ut = 
0), the length scale induced by the trap behaves as ^ 
, and the correlation function of the order parameter as 

Gix, y) = Wx)^iy))c = r'^Gixr' ,yl-''). (3) 

Finite size effects, due to a finite volume L^, can be taken 
into account by adding a further dependence on Ll~^ in 
the above scaling Ansatz [l5| . 

The value of 9 depends on the way the external con- 
fining field is coupled to the model variables. In the 
case relevant for the above-mentioned particle systems, 
the external trapping potential is coupled to the parti- 
cle density. The corresponding perturbation can be in- 
ferred from the many-body Hamiltonian [l6| in the pres- 
ence of the external potential V{x) = v^x^, i.e., Py = 
J d'^xV{x)\ip{x)\'^ , where tp{x) is the complex order- 
parameter field. The computation of the RG dimen- 
sions of the trap parameter leads to 0] = pv/{l +pi^)- 
In order to apply it to the KT transition of 2D U(l)- 
symmetric systems, we formally set — oo, correspond- 
ing to the KT exponential behavior of the correlation 
length ^ - exp(r-i/2) where r = T/T^ -1^0+, thus 
obtaining 9 — 1 for any power p of the potential. For 
comparison, we mention that 9 = 0.57327(4) at the 3D 
BEC transition in a harmonic trap Moreover, in a 
Gaussian theory perturbed by Py, since v = 1/2 we have 
9 = p/ {2+p), for any spatial dimension. It is worth men- 
tioning that analogous TSS behaviors with the same 
trap exponents, apply to the quantum T = superfluid- 
Mott transition of the BH model ([T]) at fixed integer den- 
sity, which belongs to the ((i+l)-dimensional XY univer- 
sality class . 
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which presents the same universal features of 2D sys- 
tems whose phase diagram shows a KT transition be- 
tween high-T and low-T QLRO phases. We may further 
exploit universality to investigate the effects of an in- 
homogeneity analogous to that of 2D trapped particle 
systems, where the external confining potential is gener- 
ally coupled to the particle density, which can be asso- 
ciated with an energy-density operator in a correspond- 
ing effective model. A 2D XY model with an external 
space-dependent field coupled to the energy density is 
obtained by considering a space-dependent hopping pa- 
rameter Uij, 17 1 



Hu ^ -j'^RetpiUijipj, (5) 

U,j = l + V{uj), Vir)^vPrP, (6) 

where p is an even positive integer, r.y is the distance 
from the origin of the midpoint of nearest-neighbor sites. 
We set J = 1. The inhomogcneity arising from the space 
dependence of Uij is analogous to that arising from a 
trapping potential in particle systems, such as the BH 
model ([1]). Thus, I = 1/v may be considered as the 
analog of the trap size. At large distance, since V{r) ^■ 
oo, the spin variables get effectively frozen. When p 
oo the effect of the external potential V is equivalent to 
confining a homogeneous system in a box of size L = I = 
1/v, and the TSS becomes the standard finite-size sca ling 
(FSS). At the critical KT temperature = 0.893(1) [l8[ 
of the homogeneous model (j4|) , the TSS of the model ^ 
is expected to follow the scaling Ansatz ((2]) and ((3]), with 
T] — 1/A and 9 = I as computed above by RG arguments. 

We now turn to the QLRO phase, which is the main 
issue of this paper. The homogeneous model is criti- 
cal in the whole low-T region T < Tc, where the cor- 
relation function {'tpxi'y) decays as 1/\x — yl'''"^-* with a 
T-dependent exponent r^{T): ri(T) = T/{2tt) + 0{T'^), 
increasing up to rjiTc) = 1/4 (some numerical estimates 
are reported in Refs. 0, [l^]). This critical behavior is 
controlled by a line of Gaussian fixed points, essentially 
given by the spin- wave theory Hgw = / d^x (V(^)^, which 
is the leading nontrivial term for T — > 0. We may apply 
the same spin-wave approximation to infer the value of 9 
controlling the TSS in the QLRO phase. In the spin- wave 
limit we obtain 



i^sw= / d'xl{l+vPrP){V^y^ 



(7) 



The trap exponent 9 is related to the RG dimension yy 
of parameter v, 9 = 1/yv, which can be obtained from 
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FIG. 1: Log- log plot of mo = 
T = 0.4, 0.5, 0.8, for L/i ~ 2 [H 
visible. The lines show fits to al 
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FIG. 2: G{0,x)/ml vs x/l for the model ©, at T = 0.5, 0.8, 
for several values of I, and L/Z ~ 2 [SJ. The sets of data at 
fixed T are clearly converging to a nontrivial large-Z limit 



the relation pj/t, — p + = d taking also into account 

that y(d,ip)^ = in the spin- wave theory. We eventually 
obtain 9 = 1 independently of the power p. 

Therefore, we have = 1 at very low-T, where the spin- 
wave approximation holds, and at the KT transition at 
T = Tf.. A natural scenario is that the TSS has universal 
features in the whole QLRO, with 6* = 1 for any T < Tc- 
We shall provide numerical evidence of this scenario. 

For this purpose, we present numerical results for the 
model ^ with V{r) = (r/0^ for T < ^ 0.893, ob- 
tained by Monte Carlo (MC) simulations. We center the 
trap in the middle of a square lattice (2L -I- 1) x (2L -1-1). 
The parameter L and I are chosen so that the spin vari- 
ables close to the boundaries are effectively frozen, mak- 
ing unnecessary the use of larger lattices. We use fixed 
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boundary conditions = 1- 

We consider the local magnetization at the origin and 
the two-point correlation function, defined as 

mo = {ipo), G(f,27) = (V^sVa) - (8) 



The boundary conditions = 1 breaks the U(l) sym- 
metry, thus allowing a nonzero local magnetization. Ac- 
cording to the above scaling considerations, we expect 
that their asymptotic trap-size dependence is 

mo ^ r"(^)«/2, (9) 
G{x,y)^r^^^^'g{x/l',y/l'), (10) 

where 9 = 1 and rj{T) is the T-dependent exponent of the 
homogeneous system. Since 9 = 1, a nontrivial simulta- 
neous FSS and TSS limit can be achieved by keeping L/l 
fixed, where scaling behaviors analogous to Eqs. ^ and 
pUI apply. 22] 

Fig. [1] shows data of mg for some values of T < 
Tc- A power-law trap-size dependence is clearly sup- 
ported by the data. Fits to aZ"''/^, see Fig. [1] give 
C = 0.072(1), 0.093(1), 0.179(1) respectively for T = 
0.4, 0.5, 0.8. Since according to Eci. ® C = '7^' and 9 = 1, 
these results should be compared with the available esti- 
mates of ri{T) obtained in homogeneous systems, which 
are in good agreement [2^ . 

Fig. m shows results for G{0,x). In agreement with 
Eqs. (HD and they show that G(0, x)/ml = g{x/l) in 
the large-Z limit. At small distance x ^ I, g{y) y~^^'^\ 
to recover the behavior of the homogeneous system. 

The spatial dependence of the local magnetization does 
not show a simple scaling behavior. The numerical re- 
sults appear consistent with 



(11) 



where X = x/l and Tx = T/{1 + X^), which can be 
derived using arguments based on a local-temperature 
approximation, noting that Tx may be considered as an 
effective local (space-dependent) temperature. 

We also consider a 2D XY model where the hopping 
parameter decreases moving far from the origin, i.e., re- 
placing 



U.^ = [i + Vin,)]- 



V{r) 



2 2 



(12) 



in Eq. ([5]). In this case the regions far from the origin are 
effectively in the high-T phase. The lattice system is set 
as before, but we use open boundary conditions which 
are compatible with a diverging hopping parameter at 
large distance. Thus, (ipx) = everywhere, including 
the origin. MC simulations show that the TSS is again 
characterized by the trap exponent 9 = 1. This is shown 
by Fig. [HI where we report l^^'^^G{0,x) for several values 
of I, L/l w 2 jUl, with 77(T) obtained from the data 
of Fig. [H The sets of data for T = 0.5 and T = 0.8 are 
clearly converging to a nontrivial large-Z limit, at least up 
to X EE x/l corresponding to Tx c± T/(l + X'^) = Tr^ 
0.893 {Xc « 0.89, 0.34 for T = 0.5, 0.8 respectively). [2^ 
In conclusion, we have characterized the trap-size de- 
pendence within the low-T QLRO phase and at the KT 
finite-T transition of 2D trapped systems. Using scaling 
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FIG. 3: r(^)G(0,x) vs x/;, for the U^J given in Eq. dH]), at 
T = 0.5, 0.8, for several values of I, and L/l ~ 2. 

arguments, we have argued that it is described by the 
TSS Ansatz ^ and with the trap exponent 9 = 1 
in the whole QLRO phase, up to the KT transition. This 
scenario has been supported by numerical results for clas- 
sical 2D XY models with space-dependent hopping pa- 
rameters, which give rise to inhomogeneities analogous 
to that of trapped atomic gases in actual experiments. 
These results should be useful to get evidence of QLRO 
in trapped systems, and also determine the critical pa- 
rameters at the KT transition, such as Tc in model ([5]). 
For example, the KT critical point corresponds to a TSS 
with rj ~ 1/4, while values rj < 1/A corresponds to the 
low-T QLRO phase. 
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